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Is HCN a True Tracer of Dense Molecular Gas in LIRGs and 

ULIRGs? 

J. Gracia-Carpio^, S. Garcia-Burillo^, P. Planesas^ and L. Colina^ 

ABSTRACT 

We present the results of the first HCO^ survey probing the dense molecu- 
lar gas content of a sample of 16 luminous and ultraluminous infrared galaxies 
(LIRGs and ULIRGs). Previous work, based on HCN(l-O) observations, had 
shown that LIRGs and ULIRGs posses a significantly higher fraction of dense 
molecular gas compared to normal galaxies. While the picture issued from HCO"'" 
partly confirms this result, we have discovered an intriguing correlation between 
the HCN(1-0)/HCO"'"(1-0) luminosity ratio and the IR luminosity of the galaxy 
(Lir). This trend casts doubts on the use of HCN as an unbiased quantitative 
tracer of the dense molecular gas content in LIRGs and ULIRGs. A plausible 
scenario explaining the observed trend implies that X-rays coming from an em- 
bedded AGN may play a dominant role in the chemistry of molecular gas at 
-^iR > 10^^ Lq. We discuss the implications of this result for the understanding 
of LIRGs, ULIRGs and high redshift gas-rich galaxies. 

Subject headings: galaxies: active — galaxies: ISM — galaxies: starburst — 
infrared: galaxies — ISM: molecules — radio lines: galaxies 



1. Introduction 

Since their discovery by the IRAS survey, luminous infrared galaxies (LIRGs, Lir > 
10^^ Lq) and their extreme counterparts, ultraluminous infrared galaxies (ULIRGs, Ljr > 
10^^ Lq), have been studied at virtually all wavelengths (Sanders & Mirabel 1996). The 
infrared (IR) emission of LIRGs and ULIRGs is known to stem from dust reprocessing of 
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higher frequency radiation. However the nature of the power source of the huge IR emission 
in these galaxies is still under debate (Sanders et al. 1988; Veilleux et al. 1995; Genzel et al. 
1998; Imanishi et al. 2006; Risahti et al. 2006). 

LIRGs and ULIRGs are known to possess large amounts of molecular gas as derived 
from CO(l-O) observations (Sanders et al. 1991; Solomon et al. 1997; Gracia-Carpio et 
al., in prep.). The bulk of the molecular gas reservoir of ULIRGs lies in their central kpc 
(Downes & Solomon 1998). Sanders et al. (1988) first reported that the infrared-to-CO 
luminosity ratio in ULIRGs is anomalously high compared to that of normal galaxies and 
interpreted this result as evidence of the AGN power source scenario for ULIRGs. Gao 
& Solomon (2004a, b, hereafter GS04ab) used HCN(l-O) observations to probe the dense 
molecular gas content of a sample of 65 nearby galaxies, including 25 LIRGs and 6 ULIRGs. 
Their results, showing a tight linear correlation between the IR and HGN luminosities over 
3 orders of magnitude in Lir, were interpreted in terms of star formation as being the main 
power source in ULIRGs. However, doubts have been casted on the reliability of HGN as an 
unbiased tracer of dense molecular gas in LIRGs and ULIRGs. First, X-rays may significantly 
enhance HGN abundances in enshrouded AGNs (Lepp & Dalgarno 1996; Kohno et al. 2001; 
Usero et al. 2004). Furthermore the excitation of HGN lines in LIRGs and ULIRGs might 
be affected by IR pumping through a 14 /im vibrational transition (Aalto et al. 1995). Taken 
together, the possible caveats on the use of HGN observations call for the use of alternative 
tracers of dense gas in LIRGs and ULIRGs. 

In this Letter, we present observations made with the IRAM ^ 30m telescope in the 
1-0 and 3-2 lines of HGO^ of a sample of 16 galaxies including 10 LIRGs and 6 ULIRGs. 
The new HGO+ spectra are complemented with GO (1-0) and HGN (1-0) observations either 
existing or acquired for this work. The HGO+ and HGN J=l-0 lines have both similar 
critical densities and the bulk of their emission is expected to arise from dense molecular gas 
(n(H2) > lO^cm-3). Preliminary results of this HGO+ survey indicate that the HGN/HGO+ 
luminosity ratio sharply increases with Ljr for LIRGs and ULIRGs. 

2. Observations 

The observations were carried out in November 2004 and August 2005 with the IRAM 
30m telescope at Pico Veleta (Spain). We have obtained single-pointed spectra toward 16 
objects with clear detections in all cases for the 1-0 line and with 10 detections for the 
3-2 hue of HGO"*". Our sample encompasses a LiR,-range=10^^-^-10^^-^ L©. All targets are 
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distant enough so as to recover the bulk of the emission from the molecular gas disks with the 
30m telescope beam (FWHM ~28"@89GHz and ~9"@268GHz). HCN(l-O) and C0(1- 
0) single-dish data for the sample galaxies were taken from the literature (Sanders et al. 
1991; Solomon et al. 1992, 1997; GS04ab) or were observed by us with the 30m telescope 
(Table 1). To compare the luminosity ratios derived for LIRGs and ULIRGs with those of 
normal galaxies (Lir < 10^^ Lq), we have compiled a sample of 69 objects for which HCN, 
HCO+ and CO data are available (Nguyen-Q-Rieu et al. 1992; GS04ab; this work). 

3. Caveats on the Use of HCN as Tracer of Dense Gas 

ULIRGs in GS04ab's sample have HCN-to-CO luminosity ratios up to ~10 times higher 
than measured in normal galaxies. This result is corroborated by the new HCN observations 
of the sample of LIRGs and ULIRGs conducted for this work. As shown in Figure la, the 
HCN(1-0)/CO(1-0) luminosity ratio stays fairly constant around ~0. 03-0. 04 for the normal 
galaxies of GS04ab's sample. Though with considerable scatter, the corresponding ratio 
for LIRGs and ULIRGs sharply increases, ranging from ~0.05 to ~0.30. However, if the 
HCO'''(1-0) line is adopted as an alternative tracer of the dense molecular gas, the situation 
changes significantly. As illustrated in Figure lb, we find an average HCO"'"(1-0)/CO(1-0) 
ratio close to ~0.04 in normal galaxies using the sample of Nguyen-Q-Rieu et al. (1992). 
In the case of LIRGs and ULIRGs we do derive a comparatively higher HCO^/CO ratio, 
suggestive of a correspondingly higher fraction of dense gas at higher Lir. However, the 
estimated HCO"'"/CO ratio lies between ~0.05 and ~0.13, i.e., a factor of 2-3 less than the 
corresponding range derived for the HCN/CO ratio in LIRGs and ULIRGs. Further evidence 
that the dense molecular gas fraction of LIRGs and ULIRGs does not increase dramatically 
with LiR is provided by Figure Ic which shows that the derived HCO+(3-2)/CO(l-0) ratio 
is roughly constant as a function of Ljr for Ljr > 10^^ Lq with an average value of 0.026. 
We emphasize that the 3-2 line of HCO"'" is a secure tracer of dense (and warm) molecular 
gas, as the critical density for this transition is ricrit > lO^cm"^. 

Figure 2a represents the HCN(1-0)/HCO+(1-0) luminosity ratio as a function of Lir 

for the LIRGs and ULIRGs in our sample. In spite of the small sample, we find an intriguing 
correlation between the HCN/HCO"'" luminosity ratio and Ljr. The HCN/HCO+ ratio is 
seen to increase by a factor of 5 from LIRGs to ULIRGs in our sample. We discuss below 
various scenarios that can explain why HCN seems to be over-luminous with respect to 
HCO+ in LIRGs and ULIRGs: 

Infrared Pumping of the HCN(l-O) Line — It has been argued that the excitation of 
HCN (1-0) in molecular clouds could be enhanced by infrared pumping via a 14 /xm vibra- 
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tional transition near strong mid-infrared sources (Aalto et al. 1995). LIRGs and ULIRGs 
are likely to develop this type of environment for molecular gas. To evaluate infrared pump- 
ing we should ideally know the intensity of the mid-IR field felt by the molecules in each 
source. This would require a number of assumptions to be made on the geometry of the 
problem and on the degree of dumpiness of molecular gas. Following an empirical approach, 
we have compared the ratio of HCN and HCO^ luminosities with an estimate of the excess in 
the 12 fj,m IRAS band: the ratio between the 12 //m and 100 //m IRAS density fluxes fi2/fioo 
(see Figure 2b; the 12 /im IRAS band covers the range of the HCN vibrational transition at 
14 //m). The observed lack of correlation between HCN(1-0)/HCO+(1-0) and fi2/fioo argues 
against infrared pumping of HCN(l-O) (see also discussion in GS04ab). However, caution 
is required when interpreting the fi2/fioo ratio as a measure of the MIR excess in the case 
of heavily obscured LIRGs and ULIRGs due to the high extinction in the MIR range (e.g., 
Soifer et al. 2002; Imanishi et al. 2006). The absence of correlation can not be taken as a 
firm evidence against IR pumping. On the other hand, the HCO^ line could also be affected 
by IR pumping via a 12 //m vibrational transition. Even more, the conditions for infrared 
pumping for HCN (1-0) and HCO"'"(1-0) are similar in terms of the required intensity of the 
IR field (Carroll & Goldsmith 1981), a result that argues against this scenario. 

Chemical Enhancement of HCN in Star Forming Regions — There is observational evi- 
dence that the abundance of HCN can be enhanced in the molecular gas closely associated 
with high- mass star forming regions (e.g., Blake et al. 1987). In contrast with quiescent 
molecular clouds, the measured HCN(1-0)/HCO+(1-0) luminosity ratio is high (>1.5-2) in 
galactic GMCs (Turner & Thaddeus 1977). More recently Pirogov (1999) published a survey 
of 34 bright FIR sources in the outer Galaxy identified as embedded high-mass star forming 
cores. Their results confirm that the HCN(1-0)/HCO"'"(1-0) luminosity ratio is high in these 
regions: we have estimated that the average HCN(1-0)/HCO+(1-0) luminosity ratio is ~2.3 
for their cores sample. The densest (n(H2) > 10^ cm^^) and hottest (T^ > 100 K) phase of 
the molecular star forming complexes (the 'hot core' phase) can be a privileged niche for 
HCN chemistry, favoured by the evaporation of grain mantles. If the number of hot cores 
in LIRGs and ULIRGs exceeded that in our Galaxy by at least an order of magnitude, the 
observed HCN(1-0)/HCO+(1-0) luminosity ratio could be thus accounted by a pure star 
formation scenario. In this case the reliability of HCN as a straightforward tracer of dense 
molecular gas should be equally put on hold. 

Chemical Enhancement of HCN in X-ray Dominated Regions (XDR) of AGNs — Molecular 
gas can be exposed to a strong X-ray irradiation close to the central engines of active galax- 
ies. In contrast to UV-photons, hard X-rays are efficient at penetrating huge gas column 
densities out to Ay = 100-1000 (Lepp & Dalgarno 1996; Maloney et al. 1996). There is 
founded evidence that the circumnuclear disks (CND) of some nearby Seyferts can become 
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giant X-ray dominated regions (Kohno et al. 2001; Usero et al. 2004). In the case of ULIRGs, 
it has been proposed that buried AGNs may create extended XDR rather than develop Nar- 
row Line Regions (NLR) around their central engines, making the identification of AGNs in 
the optical difficult (Imanishi et al. 2006). In particular, the influence of X-rays can enhance 
the abundance of HON relative to other tracers of the dense molecular gas such as HCO"'" 
(Lepp & Dalgarno 1996). 

Kohno et al. (2001) have proposed the use of a diagnostic tool to distinguish between 
'pure' AGNs and 'composite' starbursts+AGNs in nearby Scyfcrts, by exploring the range 
of variation of the HCN(1-0)/HCO+(1-0) and HCN(1-0)/CO(1-0) ratios. Wc have applied 
such diagnostic tool to our sample. As shown in Figure 2c, LIRGs and ULIRGs mostly 
populate the location of the 'composite' AGN family in the diagram. However, compared 
to nearby Seyferts, we find more LIRGs/ULIRGs in the transition zone between 'composite' 
and 'pure' AGNs (see Figure 3 of Kohno 2005). This is an indication that XDR chemistry 
could be at work in enshrouded AGN in these intermediate objects. According to Figure 2c, 
Mrk 231 and Mrk 273 happen to lie in the 'pure' AGN area. Of particular note, these 
two ULIRGs present firm evidence for hosting an embedded AGN in other wavelengths (see 
Table 1). Contrary to Kohno et al. (2001), who derived their ratios from interferometer 
intensity maps of the CND of nearby Seyferts, the ratios shown in Figure 2c for LIRGs 
and ULIRGs represent total disk luminosity ratios. Therefore the case for a giant XDR in 
Mrk 231 and Mrk 273 is more compelling than in nearby Seyferts as their ratios are here 
weighted by the total content of their molecular disks, including the contribution from the 
star burst and the AGN. 



4. Is HCO+ a True Tracer of Dense Gas in ULIRGs? 

Most of the theoretically founded doubts have been thus far casted on the interpretation 
of the HCN/CO ratio as to provide a quantitative estimate of the fraction of dense molecu- 
lar gas in ULIRGs. However, the intriguing trend shown in Figure 2a could be alternatively 
interpreted as due to a possible deficiency of HCO+ relative to HON in ULIRGs. In particu- 
lar, it is tempting to try to make a link between a hypothetical deficiency of HCO"'" and the 
well known deficiency of the [C II] 158 yum line relative to Lfir in ULIRGs (Luhman et al. 
2003). Observational and theoretical evidence argue against this scenario, however. First, 
the bulk of the [CII] 158 /xm and HCO+(1-0) line emissions arc not expected to arise, on 
average, from the same neutral phase of the ISM. Besides the fact that the critical densi- 
ties of the two transitions markedly differ (ucrit HCO"'"(1-0) > lOxncrit [CII]), PDR models 
predict that the bulk of the contribution to low-J HCO"*" fines should come from the inner 
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layers of the PDR where the abundance of 0+ is diminished (Sternberg & Dalgarno 1995). 
Moreover, significant non-PDR contributions to the HCO"'"(1-0) hne are very hkely. In addi- 
tion, most of the explanations brought forward thus far to solve the [C II] puzzle in ULIRGs 
invoke a quenching of the 158 fim line of C+ rather than a lower ionization degree in the 
PDR phase anyway (Luhman et al. 2003). In conclusion, there is no well founded evidence 
pointing to a physical link between a hypothetical deficiency of HCO"'" in ULIRGs and the 
well-estabhshed deficiency of [C II] 158 /xm line in these galaxies. Other possible scenarios 
invoking a deficiency of HCO+ should be explored in the future, however. 



5. Conclusions and Perspectives 

Preliminary results of the first HCO^ survey conducted in LIRGs and ULIRGs provide 
evidence that HCN may not be a true tracer of dense molecular gas in these galaxies. Differ- 
ent mechanisms can make for HCN(l-O) being over-luminous with respect to HCO+(1-0). 
The caveats on the interpretation of HCN observations highhght the need of surveys in other 
molecular species that together provide an unbiased estimate of the dense molecular gas frac- 
tion of LIRGs and ULIRGs. This question is central to disentanghng the different power 
sources of the huge infrared luminosities of these galaxies. The variation of the HC0"'"(1- 
0) /CO (1-0) ratio among the sample galaxies shows that the fraction of dense molecular gas 
of LIRGs and ULIRGs is, on average, a factor of ~2 higher than that of normal galaxies 
(Figure 16). This is considerably less than the corresponding number derived from HCN 
(i.e., a factor of ~4-5). The increase in the dense gas fraction derived from HCO"*" would 
fall short of explaining the observed Ljr for ULIRGs in the purely star formation scenario of 
GS04ab. Instead, this result suggests that the contribution to Lir from an embedded AGN 
source would amount to ~50% in the most extreme ULIRGs of our sample (Mrk 231 and 
Mrk 273). Most remarkably, the correlation found between the HCN(1-0)/HCO"'"(1-0) ratio 
and LiR for LIRGs and ULIRGs runs in parallel with the long known general tendency of 
finding more AGN signatures in ULIRGs with increasing Ljr (Veilleux et al. 1995). This 
seems to be confirmed by the derived location of Mrk 231 and Mrk 273 in the correlation plot 
of Figure 2a, but also by the position of these galaxies in the diagnostic tool of Figure 2 c. 
More recently, HCN observations of an ultraluminous quasar a,t z — 3.9 (Wagg et al. 2005) 
have shown that the HCN/CO abundance ratio could be increased due to the presence of an 
AGN that dominates the bulk of the IR emission. 

The extension of this work to galaxies with higher Lir will help to shed light on the 
relative contribution of star formation and AGN to the huge infrared luminosity of high-z 
objects. In particular, future results coming out from the extension of this survey to other 
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molecular lines will allow us to better constrain the physical and chemical status of the dense 
molecular gas in ULIRGs (Gracia-Carpio et al., in prep.). 

This work has been partially supported by the Spanish MEG and Feder funds un- 
der grant ESP2003-04957 and by SEPGT/MEG under grants AYA2003-07584 & AYA2002- 
01055. 
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Fig. 1.— Variation of the HCN(1-0)/CO(1-0) (a) and HCO+(1-0)/CO(1-0) (b) luminosity 
ratios from normal galaxies (Lir < 10^^ Lq) to LIRGs and ULIRGs (Lir > 10^^ Lq). The 
open circles (triangles) correspond to galaxies from GS04ab's (Nguyen-Q-Rieu et al. 1992) 
sample. The HCO+(3-2)/CO(l-0) ratio (c) shows little variation within the explored range 
for LIRGs and ULIRGs. Filled circles correspond to our sample in all panels. 
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Fig. 2. — a) The correlation between the HCN(1-0)/HCO"*"(1-0) luminosity ratio and Lir is 
displayed for the sample of LIRGs and ULIRGs considered in this paper (the linear regression 
fit to the data gives a correlation coefficient R = 0.84). b) The HCN(1-0)/HCO+(1-0) ratio 
shows no correlation with the infrared colour fi2/fioo- c) The location of LIRGs and ULIRGs 
of our sample in the 2D diagnostic diagram of Kohno et al. (2001). Filled symbols represent 
LIRGs and ULIRGs with secure identification of an embedded AGN. 



Table 1. Sample Properties and Observational Results 



Dl i-iR i'CO(i-o) -^-HCNCi-o) -^-HCO+Ci-o) ■'^HCO+(3-2) Nuclear Classification 



Source 


Mpc 10" Lq 


10* L' 


10* L' 


10^ L' 


10» L' 


L-band MIR Opt. X-ray Ref. 


Mrk 231 


170 


29.3 


74 


20 


8.2 


2.5 


AGN 


AGN 


Syl 


AGN 1 


IRAS 17208-0014 


173 


21.9 


131 


16 


11 


3.7 


SB 


SB 


HII 


SB 1 


IRAS 12112 + 0305 


297 


17.6 


101 


13 


7.6 




SB 


SB 


LIN 


SB 1 


IRAS 23365 + 3604 


261 


13.2 


69 


6.6 


4.1 






SB 


LIN 


SB 2,3,4 


Arp 220 


73 


13.0 


70 


8.3 


4.8 


1.1 


SB 


SB 


LIN 


AGN? 1 


Mrk 273 


152 


12.5 


55 


13 


5.7 


2.8 


AGN 


AGN 


Sy2 


AGN 1 


VII Zw 31 


220 


8.2 


114 


8.4 


6.9 








HII 


5 


NGC 6240 


98 


6.1 


82 


13 


9.6 


2.0 


AGN 


SB 


LIN 


AGN 1 


Arp 55 


161 


4.5 


127 


5.1 


4.2 








HII 


2 


NGC 695 


131 


4.2 


48 


1.7 


2.6 






SB 


HII 


2,6 


Arp 193 


94 


3.9 


38 


1.8 


3.5 


1.3 




SB 


LIN 


2,6 


NGC 7469 


66 


3.9 


34 


2.0 


2.9 


0.35 




SB 


Syl 


AGN 2,7 


Arp 299 A 


42 


3.2 


13 


0.89 


1.7 


0.37 




SB 


HII 


AGN? 5,7,8 


Arp 299 B+C 


42 


2.5 


14 


0.54 


0.84 


0.11 




SB 


HII 


AGN 5,7,8 


Mrk 331 


72 


2.6 


42 


1.9 


2.4 


0.34 




SB 


HII 


••• 2,6 


NGC 7771 


57 


2.1 


43 


3.2 


3.0 


< 0.70 




SB 


HII 


••• 2,6 



Note. — The luminosity distance, Dj,, was calculated assuming Ho = 75 Kms-l Mpc"! and go = 0.5. The 
luminosities were computed according to Sanders & Mirabel (1996) and Solomon et al. (1997) formulas. L' = K 
kms^^ pc^, in Tmb scale. Single-dish data taken from the literature arc marked in boldface. References for the 
nuclear classifications: (1) Risaliti et al. 2006, and references therein; (2) Veilleux et al. 1995; (3) Rigopoulou et al. 
1999; (4) Teng et al. 2005; (5) NED; (6) Lu et al. 2003; (7) Genzel et al. 1998; (8) Ballo et al. 2004. 



